Effective thermal dynamics following a quantum quench in a spin chain 
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We study the nonequilibrium dynamics of the Quantum Ising Model following an abrupt quench 
of the transverse field. We focus on the onsite autocorrelation function of the order parameter, and 
extract the phase coherence time Tq from its asymptotic behavior. We show that the initial state 
determines Tq only through an effective temperature set by its energy and the final Hamiltonian. 
Moreover, we observe that the dependence of Tq on the effective temperature fairly agrees with that 
obtained in thermal equilibrium as a function of the equilibrium temperature. 
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A recent series of beautiful experiments with cold 
atomic gases [1 H Q have triggered a great deal of inter- 
est in some fundamental aspects of the non-equilibrium 
dynamics of correlated quantum systems. The pecu- 
liarity of the dynamics of cold atoms is its phase co- 
herence on long time scales. This was clearly demon- 
strated by the cycles of collapse and revival of the order 
parameter observed in Ref. [2|. The interplay between 
phase coherence, strong interactions, and low dimension- 
ality may result in surprising effects: an example is the 
lack of thermalization recently observed in quasi-one di- 
mensional condensates [3[. The attribution of this phe- 
nomenon to the closeness of these systems to integrabil- 
ity spurred an intense discussion on the general relation 
between quantum integrability and thermalization in the 
long-time dynamics of strongly correlated quantum sys- 
tems 0,1,1101©, El El El El El El- 

The simplest nonequilibrium process to be considered 
in order to study the long-time dynamics of a quan- 
tum system is the quantum quench: an abrupt change 
in time of one of the system parameters or of its bound- 
ary conditions. Recent studies of stron gly correlated 
models @, |, @, 0, 0, ©, El EH, El ElM El have 
demonstrated that the behavior of integrable and non- 
integrable systems can be quite different. Thermaliza- 
tion can be observed, under specific circumstances, in 
nonintegrable systems d d E3] : asymptotic values of 
significant observables, such as the momentum distri- 
bution function, do not depend on the details of the 
initial state, but only on its energy 0. On the other 
hand, for integrable systems thermalization does not oc- 
cur @, H 0, El El El El EE|: a larger amount of in- 
formation on the initial state seems necessary to predict 
the asymptotic state. It has been conjectured that this 
information consists of the expectation value of a set of 
constants of motion fixing in the Lagrange multipliers of a 
generalized Gibbs ensemble Ml . For a special quench in a 
ID Bose- Hubbard model II ll and for integrable systems 
with free quasiparticles [l2|, the local reduced density 
matrix was indeed proven to asymptotically tend to such 
generalized ensemble. Moreover, the generalized Gibbs 



ensemble was shown to correctly predict the asymptotic 
momentum distribution functions for a variety of mod- 
els and quenches [1 0, E! 14 1 . However, it should be 
pointed out that neglection of correlations of the occu- 
pation of different quasi-particle modes leads to incorrect 
predictions for the noise and higher order correlators [H[ . 

In this Letter, instead of focusing on the asymptotics 
of observables, we take a different perspective, and study 
the dependence on the initial state |V>o) of the intrinsic 
time-scale of the dynamics after the quench. We do this 
by considering the Quantum Ising chain, a prototypical 
example of exactly solvable model with a quantum phase 
transition j7i§] . We study the autocorrelation function of 
the order parameter after a quench of the transverse field, 
extracting the phase coherence time Tq from its asymp- 
totic exponential decay. We will show that, regardless of 
the integrability of the model, the only information on 
the quench needed to predict Tq is the final gap A and 
an effective temperature T e ff, determined by the energy 
of the initial state after the quench; moreover, the depen- 
dence of Tq on T c g is in very good agreement with that 
obtained, at equilibrium, for the same quantity t!£ as a 
function of the equilibrium temperature T . 

The sharp contrast between the asymptotics of observ- 
ables like the transverse magnetization, determined by 
the entire set of constants of motion, and the phase co- 
herence time Tq, depending just on T e g, has its deep 
roots in the physics of the Quantum Ising chain [1 €31 ] - 
This model, which can be diagonalized in the continuum 
limit in terms of Majorana fermion quasi-particles (see, 
e -g- [I3|)i possesses two sectors of operators E3 : a lo- 
cal sector with respect to the quasiparticles, where the 
S'-matrix is simply S = 1 and the model is equivalent to 
a free theory, and a non-local sector, where S = — 1 and 
the model describes an interacting theory. In this respect, 
the non-local sector can be used as low-energy theory of a 
more general class of models, not necessarily integrable, 
belonging to the Ising universality class (e.g., the non- 
integrable <I> 4 Landau-Ginzburg model [19(). While the 
transverse magnetization belongs to the local sector, im- 
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plying sensitivity of its asymptotics to integrability, the 
order parameter belongs to the nonlocal sector, making 
Tq representative of the Ising universality class. 

We start by considering a spin-1/2 Quantum Ising 
chain in a transverse magnetic field T with periodic 
boundary conditions: 



IP 



w(r) = -j£[<#$ + i + r*S 
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where er" (a = x, y, z) are spin operators, J is the inter- 
action strength. Hereafter, unless explicitly written, we 
set J = 1. This system has a quantum critical point at 
r c = 1 separating two mutually dual gapped phases, a 
quantum paramagnetic one (T > 1) and a ferromagnetic 
one (r < 1), with energy gap A = 2|1 — T|. At equilib- 
rium, the presence of a critical point dramatically influ- 
ences the temperature dependence of the basic time-scale 
characterizin g th e system's dynamics: the phase coher- 
ence time [16( . The latter is usually extracted from the 
asymptotics of the on-site spin autocorrelation function 
Px x (t) = (o|(t)crj(0)), which decays to zero exponen- 
tially, [HElI p x T x {t) ~ e-^ T T, at any finite temperature 
T > 0, both at criticality (A = 0), and in the off-critical 
region (T < A). At criticality [2(|, for T < J one finds 
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FIG. 1: (color online). Time dependence of \pQ(t)\ for a 
quench to a final ferromagnet V — 0.5. Different curves, ob- 
tained numerically for a finite chain of L — 600 sites, corre- 
spond to different initial IYs. Inset: \pQ ! {t)\y/t for Y = 1.25. 



Most importantly, we find that Tq — t£ =t ff , both for 
quenches at criticality and away from it. 

To calculate Pq{t) we exploit the complete integrabil- 



while Tji is exponentially larger [21|, |22| in the ity f tne Ising chain g ^ . Here we sketch the essential 



off-critical region with T <C A: ~ ^e A / T . 

Consider now a quantum quench, which consists of 
preparing the system in the ground state correspond- 
ing to a transverse field Tq, |-0o) = ^(To)), and then 
abruptly quenching it, at t = 0, to some T ^ T - For 
t > 0, the state evolves unitarily under H(T), accord- 
ing to \ijj t ) = exp[— iH(T)t] |V>(r )). We define the zero- 
temperature on-site autocorrelation function describing 
the spin dynamics after the quench: 



steps 25]: first, one represents spins in terms of Jordan- 
Wigner fermions q = cr ; — exp (iix YljZ. \ (T^crj). Since the 
ground state has always an even number of fermions, one 
can focus on the even c-fermionic Hilbert space sector. 
Switching to momentum representation, the Hamilto- 
nian is diagonalized with a Bogoliubov rotation: 7i(r) = 
Efc>o e fc (^hk + l- k l-k - l), where j k are fcrmionic 
quasi-particle operators, 



2 Vr 2 - 2r cos k + 1 



their dispersion, and k = ± 



7t(2k+1) 



with n = 0, 



\t) = {i>(T )\ eWWofe-^ |V(r )) 



1. 



(2) 



Before entering into details, we summarize the results ob- 
tained by analyzing the asymptotics of pQ X (t): it always 

drops exponentially to zero (see Fig.[T]), p^it) ~ e - *^ 7 "^, 
as in the finite-temperature equilibrium case, consistent 
with what was obtained in Ref. [6( for critical quenches. 
This allows us to extract a time-scale Tq characterizing 
the dynamics after the quench. This phase-coherence 
time depends in principle on the initial state ^(To)) and 
the final Hamiltonian Ti.(T). However, and this is the 
main result of this Letter, all the information needed to 
characterize Tq is encoded in two variables only: the final 
gap A(r), and an effective temperature T c g. The latter 
is obtained by comparing the energy associated to the 
initial state with respect to the Hamiltonian after the 
quench to the average energy of a fictitious thermal state 
at temperature T c g in an effective canonical ensemble: 

e(t ) = (V(r )| w(T) |v(r )) = (w(r)>T eff . (3) 



The second step consists in describing the dynamics af- 
ter a quench. This can be easily done in the Heisen- 
berg picture (26|, by solving the closed set of equations 
of motion for the c-fermions in momentum space, with 
the initial conditions associated to the quench. Finally, 
/Og x (f) is computed using a trick developed in Ref. [25j |. 
The operator <Tj(t)aJ(Q) connects states with different 
c-fermion parity, and it cannot be simply evaluated us- 
ing Jordan- Wigner fermions in the even Hamiltonian sec- 
tor. This problem can be circumvented by considering 
a four-spin correlation function on a chain of length L, 
C x {t-L) = (of (t) af(0) o-l +1 (t) <r| +1 (0)). This correla- 
tor conserves the c-fermion parity, and can be written as 
the square root of a Pfaffian 25j, using the techniques 
of Ref. [23|. One finally recovers pffit) using the clus- 
ter property [pcfit)] 2 = hniL^oo C x (t;L), by taking the 
square root of C x (t) in the limit of large number of spins. 

As anticipated above, the zero-temperature quench au- 
tocorrelation /0q x (<) always relaxes exponentially to zero 
(see Fig.[I|, irrespective of the initial state \ip(To)) and of 
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FIG. 2: (color online). Phase coherence time Tq as a function 
of the initial transverse field To. The different curves refer 
to a ferromagnetic (r = 0.5, red circles), a critical (r = I, 
blue triangles) and a paramagnetic (F = 1.25, green squares) 
quench dynamics. Inset: effective temperature T e g vs. Fo, as 
extracted from Eq. (3), for the same values of T. 




FIG. 3: (color online). Phase coherence time Tq as a function 
of T e ff. Different symbols are for various values of F. Empty 
symbols correspond to an initial field To < T, while filled ones 
are for To > F. Straight lines denote the finite-temperature 
equilibrium values of t^_ t . The inset shows the finite-size 
scaling of Tq at criticality. Data in the main panel are for 
L = 600 (except for T = 1 and T eS < 0.07, where L = 1000). 



the final transverse field T f Tq: p Q x (t) ~ e~ t/r Q . This 
is in sharp contrast with the zero-temperature equilib- 
rium autocorrelation p^ x =0 (t), which decays as + C/t 
for T < 1, M x = (1 — T 2 ) 1 / 8 being the spontaneous 
magnetization 25]. Quenching to the paramagnetic 
side, the exponential drop is superimposed to an oscil- 
latory power-law decay. This is once again reminiscent 
of finite-temperature equilibrium case, where Pjf(t) ~ 
K{t')e~ t / T T , K(t) being the quantum zero-temperature 
correlator [2l|, which oscillates and decays as i -1 / 2 . In- 
deed, for a quench to T > 1, rescaling PQ X (t) with the 
zero-temperature factor i -1 / 2 , we recover exponential re- 
laxation (inset of Fig. [IJ . 

We now analyze the coherence time as a function of dif- 
ferent initial and final conditions. In Fig. [21 Tq is plotted 
as a function of initial To, for several final T's. A dramatic 
increase of Tq as T — > T is observed: the less the system 
goes out-of-equilibrium, the slower is the relaxation. If 
r = r, the exponential decay turns into a power-law, as 
in the zero-temperature equilibrium case, and Tq — > oo. 
The analogy with the equilibrium finite-temperature be- 
havior, where the lower is the temperature T the longer is 
Tj., is evident. It is therefore tempting to relate the two 
cases, by introducing an effective temperature T e g for 
the out-of-equilibrium system. We define T e g by com- 
paring the energy of the initial state \ip(To)) with that 
of a fictitious thermal state as in Eq. ([3]), with a thermal 
energy (H(T)) Totf = E/=>o e fc K( T cff) + n_ fc (T cff ) - 1) 
determined by an effective canonical ensemble Fermi dis- 
tribution function nk(T c g) = (l + e efc / Toff ) _1 of the quasi- 
particles 7&. A plot of T e g as a function of To, for different 
values of T, is shown in the inset of Fig. [21 Notice that, 



for each T, there are two values of Tq for which T c g is the 
same, one for Tq < T and one for Tq > T. 

The effective temperature T e g , together with the quasi- 
particle gap A at the final T, univocally determines the 
phase coherence time Tq. Numerical evidence is shown 
in Fig. [31 points with equal A(T) have the same Tq if 
effective temperatures are the same, even if Tq and T are 
different. Since the system is closed, it would be tempting 
to substitute T e ff with just the initial energy E(Tq). This 
is not always possible. Indeed, two quenches having dif- 
ferent initial energy E(Tq), but equal final gap A(r) and 
equal T e g, will exhibit the same Tq (see data for T = 1.25 
and 0.75 in Fig. [3]). For example, the two quenches 
1.21 -► 1.25 and 0.715 -► 0.75 have equal T cS ~ 0.113 
within 0.7%, (corresponding to Tq ~ 1170 within 0.3%), 
although their energies differ by 25%. However, if the 
final r is fixed, the canonical and microcanonical effec- 
tive ensembles are equivalent, i.e., E(Tq) and T e g can 
be interchanged. We also notice that, with a good ac- 
curacy, Tq is still given by the equilibrium expressions 
at temperature T e g, i.e., Tq ~ t t=t s ou ^ °f criticality. 
In the critical tendency to follow the equilibrium 

expressions is observed only at low temperatures, using a 
finite-size scaling (see inset of Fig. [3]). This is necessary 
because, while the comparison is expected to work better 
at low T e ff, the consequent importance of the long wave- 
length modes makes finite-size effects more pronounced. 

The time dependence of correlators is influenced by 
quasiparticle propagation, as exemplarily shown in a 
series of studies 0, 0, Q focusing on the asymp- 



4 




FIG. 4: (color online). Phase coherence time at F = 0.5 
computed according to the quasiparticle distribution fk (tq, 
blue squares) and to the effective quasiparticle distribution 
Hfc(Teff) [jf,, red circles). 



temperatures, in agreement with the cusp singularity of 
Fig.H 

In conclusion, we studied the phase coherence time Tq 
after an abrupt quench of the transverse field in a Quan- 
tum Ising chain. We have shown that, irrespective of 
the integrability of the model, Tq depends only on the 
quasiparticle gap A and the quasiparticle effective tem- 
perature T G ff, and provided numerical evidence of the fact 
that the dependence of Tq on T c g is close to the one ob- 
tained at equilibrium as a function of the equilibrium 
temperature. The realization of the dynamics of artifi- 
cial quantum spin chains using bosonic atoms in optical 
lattices [13] represents a concrete possibility to check our 
theoretical scenario with the available experimental tools. 

We thank N. Andrei, E. Altaian, R. Fazio, V. 
Oganesyan, V. Kravtsov, and A. Polkovnikov for discus- 
sions. G.M. acknowledges hospitality at Galileo Galilei 
Institute where this work was completed, and grants IN- 
STANS (from ESF) and 2007JHLPEZ (from MIUR). 



totics of equal time correlators at different space-points, 
p xx (r,t) = (<rf +r (t)af(t)). A similar picture, elucidating 
some of the results so far obtained and applicable in the 
off-critical re gim es (T e g <C A), can be formulated in the 
spirit of Ref. 2l|. When the transverse field is quenched, 
the initial condition for the time evolution consists of 
a state with a finite density of quasiparticles (relative 
to H(T)), characterized by a dispersion e£. For T < 1, 
these can be seen as kinks propagating with momentum 
k and velocity Vk = de\/dk. The correlator Pqit) is 
determined by the number of kinks passing through a 
sin gle site in the interval [0, t\. A combinatorial analy- 
sis [21| , together with an average over momenta, leads to 
plf(t) ~ exp[-t/f|], with (f*)- 1 = f J2k M fk, where 
fk = {ipo\ lllk \ipo) is the occupation of quasiparticle 
modes. Analogous arguments can be presented for T > 1, 
giving plf{t) ~ K(t)exp(-t/T$), where K(t) « t- 1 ' 2 is 
the equilibrium zero-temperature correlator. Note that 
here integrability is not necessary: the same reasoning 
applies to the low-energy kink of non-integrable <E> 4 the- 
ory that belongs to the same Ising universality class. 
While this picture explains the exponential decay to- 
wards zero, it is also important to notice the following: 
although the quasiparticle distribution function fk de- 
termined by the initial state and the effective thermal 
quasiparticle distribution function nk{T c g) are typically 
very different, phase coherence times Tq computed from 
fk or according to (t^)" 1 = \ J2k \ v k\ n k(T cS ) are very 
close (see Fig. SJ. In other words, we explicitly checked 
that, fixing an effective temperature through Eq. (|3|) 
or using Tq — f^, leads to qualitatively and quantita- 
tively similar (to a few percent accuracy) results, for 
A T ffj moreover, by imposing fg = we find 

f eff ~ 2A/ ln[A/(r- r ) 2 ] far from criticality and at low 
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